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Abstract 

We experimentally investigate the dynamics of water cooled from 
below at 0°C and heated from above. Taking advantage of the un¬ 
usual property that water’s density maximum is at about 4°C, this 
set-up allows us to simulate in the laboratory a turbulent convective 
layer adjacent to a stably stratified layer, which is representative of 
atmospheric and stellar conditions. High precision temperature and 
velocity measurements are described, with a special focus on the con- 
vectively excited internal waves propagating in the stratified zone. 
Most of the convective energy is at low frequency, and corresponding 
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waves are localized to the vicinity of the interface. However, we show 
that some energy radiates far from the interface, carried by shorter 
horizontal wavelength, higher frequency waves. Our data suggest that 
the internal wave field is passively excited by the convective fluctua¬ 
tions, and the wave propagation is correctly described by the dissipa¬ 
tive linear wave theory. 


Keywords: internal waves, stratified flows, Rayleigh-Benard convection 


1 Introduction 


In many geophysical and astrophysical systems, a turbulent convective fluid 
layer is separated from a stably stratified layer by a relatively sharp but de¬ 
formable interface. Examples include convective and radiative zones in stars, 
the atmospheric convective layer and overlying stratosphere, the Earth’s 
oceans, and possibly the Earth’s outer core according to the latest esti¬ 
mates of its thermal conductivity (Hirose et ah, 2013). While motions in 


the stratified layers are often neglected, these layers actually support oscil¬ 
latory motions called internal waves, which can be excited by the adjacent 
convection. The remote observation of these internal waves can be used to 


probe the interior of stars via asteroseismology (e.g. Garcia et ah, 2007). In 


addition, internal waves transport energy and momentum (e.g. Bretherton 


1969; Zahn et ah, 1997; Rogers et ah, 2013). They interact non-linearly and 


produce mean motions, responsible for, e.g., the Quasi-Biennial Oscillation 


of the equatorial zonal wind in the Earth’s tropical stratosphere (Plumb and 


McEwan, 1978), and possibly explaining the observed misalignment between 


the orbits of extrasolar planets and their host stars (Rogers et ah, 2012). 


Internal waves can also break, mixing the ambient fluid. This could have pri¬ 


mordial consequences for stellar internal composition evolution (Charbonnel 


and Talon, 2005) and result in mass loss during the late stages of massive 


stars evolution (Quataert and Shiode, 2012). Internal waves are thus essential 


for accurate models of global climate and stellar dynamics. 

Global integrated models require length scales, time scales and velocity 
scales spanning many orders of magnitude to simultaneously resolve motions 
in turbulent convective and stratified zones, and to understand the details 
of the highly non-linear couplings between convective turbulence, waves and 
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meridional circulation. This is clearly very challenging both analytically and 
numerically. Recent studies have estimated analytically the flux of internal 


waves (e.g. 

Lecoanet and Quataert 

, 201f 

)) and their influence on the angular 

velocity profile (e.g. 

Alvan et ah, 

2013) 

, using models of convective turbu- 


lence and assumptions about the excitation mechanism. The relevance of 
these models remain to be validated. In particular, in both atmospheric (e.g., 


Ansong and Sutherland, 2010) and astrophysics (e.g., Goldreich and Kumar 


1990) communities, the question remains: does wave excitation mainly come 


from interface deflections or from bulk stresses within the convective zone? 
The exponential increase in computing resources has recently allowed for the 
first three-dimensional nonlinear numerical simulations of convective turbu¬ 


lence and internal waves simultaneously (Brun et ah, 2011; Alvan et al. 


2014), but only at moderate turbulent intensity. Clearly, it would be com¬ 


putationally much more efficient to simulate the convective turbulence and 


the internal waves separately using dedicated numerical schemes (Belkacem 


et ah, 2009). Crucially, this would necessitate a comprehensive description 
of the coupling between the two zones. In this context, quantitative data 
from a well-controlled experiment, in a simplified but physically motivated 
configuration, is of great interest to validate analytic theories and numerical 
simulations. 

Several experimental studies of internal wave generation by convective 
turbulence have been published since the 60’s. Most studies focused on the 
transient case where an initially thermally stratified layer of water is suddenly 


heated from below (e.g. Deardorff et al., 1969): a convective mixed zone then 


forms at the bottom of the tank and progressively invades the whole depth 
of the fluid. This configuration was in particular investigated by [Michaeli^ 
et al. (2002) using modern techniques of flow measurements. Their focus was 


mostly on the evolution of the interface position and on the fluctuations of 
the heat flux. Furthermore, this work provides the first, and—to the best 
of our knowledge—the only experimental velocity maps of internal waves 
excited by convection, which were measured using the Correlation Image Ve- 
locimetry (CIV) technique. They observed intermittence in the internal wave 
motions, dominated by two types of flows: first, short wavelength, high fre¬ 
quency waves with velocity vectors with inclinations oiO = 55° ± 5° from the 
horizontal; and later, long wavelength, low frequency waves corresponding to 
almost horizontal motions. Using the dispersion relation for internal waves, 


uj = Vsind, 


( 1 ) 
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where lo is the wave frequency, N the buoyancy frequency and 9 the angle 
between the wave vector and the vertical (equivalently, the angle between the 


isophase lines or the velocity vectors direction and the horizontal), Michaelian 


et al. (2002) interpreted these results as follows. High frequency waves, whose 


propagation angle is surprisingly constant over the range of stratifications 
studied, are excited by intermittent plumes which emanate from along the 
bottom boundary, especially early in the experiment. Later, a large-scale 
circulation develops, which excites low frequency waves by organizing the 
previously random plume motions. 

Other experimental studies used salted water as a working fluid, focusing 
on the effects of an isolated buoyant element impacting a stratified layer. 
For instance, McLaren et al. (1973) and Cerasoli (1978) investigated the 
nature and energy transfer of internal gravity waves generated by a thermal 
(i.e. the instantaneous release of a fixed volume of buoyant fluid) using a 
particle tracking method and local conductivity measurements plus dye line 
visualization. More recently, Ansong and Sutherland (2010) experimentally 
studied the generation of internal gravity waves by a turbulent plume (i.e., 
the continuous release of buoyant fluid) impinging upon the interface between 
a uniform density layer of fluid and a linearly stratified layer. Using non- 
intrusive axisymmetric synthetic Schlieren measurements, they determined 
the fraction of the energy flux carried by the waves. They also noticed that 
the dominant wave frequency lies in a narrow range relative to the buoyancy 
frequency, with 9 e [27° — 58°] and a mean value of around 45°. Following 
Dohan and Sutherland] (2003), they hypothesized that turbulent fluctuations 


interact resonantly with the excited internal waves to optimize the vertical 
transport of horizontal momentum, which for a fixed displacement amplitude 
takes place at 9 = 45°. 

A particularly original experiment was performed by Townsend (1964), 
using the unusual property that water’s density maximum is at about 4°C. 
Consequently, in a tank filled with water with an horizontal bottom surface 
at 0°C and a hotter upper surface, the density stratification is unstable at 
temperatures below 4°C and stable at temperatures above. In a sufficiently 
tall tank, it is possible to produce a turbulent convective layer adjacent to 
a stably stratified layer in a coherent self-organizing system, which is rep¬ 
resentative of atmospheric and stellar conditions. Townsend (1964) used 


temperature measurements and dye visualization to detect internal waves in 
the stratified layer, but only found low frequency waves localised close to 
the interface. It is unclear whether there were waves propagating away from 
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the interface along a “favorite angle”, as described in the closely related con¬ 


figurations studied by Michaelian et al. (2002) and Ansong and Sutherland 

( |Mol ). 

In the present paper, we re-investigate this “4°C experiment” first in¬ 
troduced by Townsend (1964) using high precision local temperature mea¬ 
surements as well as global non-intrusive velocity measurements by particle 
image velocimetry (PIV), both in the convective and stratified zones. Our 
work builds on the preliminary results presented in the conference proceed¬ 
ing by Perrard et al. (2013). Our main objectives are to better describe the 
characteristics of the excited wave field and to provide quantitative data for 
validating analytical and numerical models. The paper is organized as fol¬ 
lows. In section 2, we present the experimental set-up and the diagnostic 
tools. Section 3 is devoted to the statistics of the temperature fluctuations. 
The main results from PIV measurements are given in section 4. Finally, 
conclusions and open questions are addressed briefly in section 5. 


2 Experimental set-up 


Our experimental set-up is presented in figure [^a). It consists of a rect¬ 
angular tank whose lateral boundaries are made of vacuum glass along the 
large sides and 3 cm thick acrylic along the short sides. The internal di¬ 
mensions are 18 cm x 4 cm in width and 34 cm in height. Given the aspect 
ratios of our tank, no significative motion takes place in the short direction. 
However, based on the velocity measurements presented in section the 
typical Reynolds number of our flow is about 110 in the convective zone and 
1.5 in the stratified zone: the flow is thus influenced by viscous boundary 
layers from the side walls, especially in the stratified zone. The lower and up¬ 
per boundaries are copper plates, whose temperatures are controlled within 
±0.05°C by two circulating thermostated baths, with a fixed temperature 
equal to 0°C and 22° C respectively. To increase the thermal insulation of 
the tank, foam plates cover all boundaries, except when side view visualiza¬ 
tion is required, for instance during PIV measurements: foam plates from 
the two large sides are then removed, hence increasing heat losses. With 
all foam plates in place, the global external convective heat transfer coeffi¬ 
cient of our set-up is about 2.9 W m~^ which would correspond to 7 cm 
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Figure 1: (a) Picture of the experimental set-up, showing the tank, the two 
thermistor probes mounted on a vertical linear stage, and the green laser 
light sheet for PIV measurements. In this picture, the front and bottom foam 
insulating plates and the top copper plate have been removed for visualiza¬ 
tion. (6) Example of a time-averaged vertical temperature profile measured 
several hours after initializing an experiment. The black (red) dot shows 
schematically the location where the longterm temperature fluctuations are 
measured in the convective zone (stratified zone); see an example of temper¬ 
ature measurements in figure (c) Schematic of the expected flow with a 
large convective cell plus turbulent fluctuations in the lower convective zone 
and internal waves propagating from the interface into the upper stratified 
zone. 
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thick acrylic walls. At the beginning of each experiment, the baths are first 
allowed to reach their assigned temperature. The tank is then filled with 
water at about 4°C up to a given depth, chosen according to estimates of 
heat flux equilibrium or according to previous experiments. The upper part 
of the tank is filled with water with a linear temperature profile from 4°C 
to 22°C, using the double-bucket technique frequently used in salted-water 
experiments (see, e.g., Osterf 1965). Note that this filling process allows us to 
avoid a very long transient (typically 1 day) before reaching a quasi-steady 
state set by thermal diffusion in the upper stratified zone. The system is 
typically allowed to evolve for 1 hour before measurements are taken. 

We have performed two types of measurements. First, two high-precision 
Platinium resistance sensors with diameter 1 mm and acquisition frequency 
2 Hz give access to temperature fluctuations at chosen locations with a preci¬ 
sion of ±0.001°C (more than one order of magnitude better than the prelim¬ 
inary measurements presented in Perrard et ah, 2013). These two probes are 
mounted on a vertical linear stage with a vertical spacing of 12 cm, positioned 
in the middle of the tank in the horizontal directions. Vertical time-averaged 
temperature profiles, such as the one shown in figure [^6), are determined by 
taking measurements at regularly spaced locations. The measured profiles are 
characteristic of a two-layer system, with a monotonically increasing temper¬ 
ature in the upper stratified zone, and a fully mixed lower convective zone of 
constant temperature at about 3.4°C in the horizontal middle of the tank, i.e. 
where rising cold structures are expected (see figure [^) . The observed depth 
of the convective zone is hconv = 20 cm for the well-insulated temperature 
measurements: this corresponds to a typical Rayleigh number Ra of about 
5 X 10^. For the less-insulated experiments involving velocity measurements, 
hconv = 9 cm, corresponding to a smaller Rayleigh number Ra = 5 x 10®. 
In both cases however, Ra is above the expected transition for turbulent be¬ 
havior (Krishnamurti and Howard, 1981). The temperature gradient in the 
upper stratified zone is not rigorously constant because of side heat losses, 
despite efforts to minimize them. Combined with the peculiar equation of 
state of water, this leads to a non-constant buoyancy frequency in the strat¬ 


ified zone, N (see also Townsend, 1964; Perrard et ah, 2013; Lecoanet et al. 


2015). However, to simplify our analysis, we approximate N to be constant 


and equal to its mean value in the 6.5 cm region extending above the 8°C 
interface, where we have measured propagating waves in our set-up. This 
value is V ~ 0.06 Hz when using full insulation, and N ~ 0.07 Hz otherwise. 
Once these mean properties of the flow are established, long measurements 
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of temperature fluctuations with time at selected locations are recorded over 
3 days, and analyzed in detail in section 3. 

In addition to temperature measurements, we also perform velocity mea¬ 
surements using particle image velocimetry (PIV). The water is then seeded 
with spherical reflective particles of typical diameter 10 jim. The tank is 
illuminated from one of the short sides using a 40 mW green laser vertical 
light sheet, and a camera records the particles displacements from one of the 
long sides at 4 fps with a resolution of 1920 x 1080 pixels for the convective 
zone, and 0.5 fps with a resolution of 3680 x 2456 pixels for the stratified 
The captured movies of typical duration 10 min are PIV-processed 


zone. 


using the freely available software DPIVSoft (Meunier and Leweke, 2003). 


The velocity fields are typically resolved into 32 x 32 pixels boxes with 50% 
of overlapping, leading to a spatial resolution of about 1.3 mm in the strat¬ 
ified zone. Figure [^c) shows a sketch of the expected velocity profile, with 
a large-scale circulation plus fluctuations in the convective zone and internal 
waves in the stratified zone, as will be demonstrated in section 4. 


3 Analysis of the temperature signals 


Figure shows temperature data recorded over 3 days at fixed locations 
in the convective and stratified zones, at respective heights of 10 cm and 
22 cm from the bottom of the tank. Despite our efforts to minimize lateral 
heat losses, the temperature in the stratified zone is still very sensitive to 
the room temperature, whose variations of about 3°C between day and night 
explain the secular variations observed in figure [^a). Figures [^6) and (c) 
show the temperature fluctuations over one hour after removing the secular 
variations, in the stratified and convective zones, respectively. Note the dif¬ 
ferent temperature scales, with convective temperature fluctuations in panel 
(c) being more than one order of magnitude larger than internal wave tem¬ 
perature fluctuations in panel (6), even though the stratified probe is located 
only ~ 2 cm above the interface. Temperature fluctuations in the convec¬ 
tive zone shown in figure [^c) are clearly skewed, with sharp negative peaks 
corresponding to cold plumes rising in the middle of the tank. On the con¬ 
trary, temperature fluctuations in the stratified zone shown in figure |^6) are 
more symmetric, as expected from temperature fluctuations associated with 
propagating internal waves. 
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(a) full signal 
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Figure 2: (a) Temperature measurements during 3 days at fixed locations 
in the convective zone (lower black signal) and in the stratihed zone (upper 
red signal) as a function of time in hours, {b) and (c) show a zoom of the 
corresponding fluctuations between 13 and 14 hours—shown by the gray 
rectangle in (a)—once secular trends have been removed by a polynomial fit. 
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This conclusion is confirmed in figure [^a), which shows the probability 
density function of each of the measured signals normalised by its maxi¬ 
mum value. To account for the longterm variations shown in figure |^a), the 
raw temperature measurements are cut into 4-hours long signals, from which 
secular variations are removed by a polynomial fit of degree 10 in the strati¬ 
fied zone and degree 1 in the convective zone. Each result shown here then 
corresponds to the average probability density function over the 18 signals 
obtained from the 3-days measurements. The statistics of the temperature 
fiuctuations in the convective region is comprised of two superimposed Gaus- 
sians: one from the convective turbulence, and the other from the ascending 
cold plumes. On the contrary, the statistics of the temperature fluctuations 
in the stratified zone shows an exponential distribution, which is due to the 
intermittent generation of internal waves. 


We also computed the power spectral densities (PSD) of the same signals, 
which are shown in figure |^6) . The upper black convective signal is typical 
of turbulent Rayleigh-Benard convection, with a decreasing negative slope 
towards high frequency. The peak at 2 x 10“^ Hz corresponds to the 
overturning timescale, and is in good agreement with the observed period 
of temperature fluctuations (see figure [^). We show for comparison the 
power laws that we estimated from the data published in Wu et al. (1990), 
obtained in a classical Rayleigh-Benard experiment at Rayleigh numbers 7 x 
10® - 1.1 X 10®. The good agreement with our measurements suggests that 
the presence of the stratified layer does not significantly modify the behavior 
of the convective layer. 

As expected, the power spectral density in the stratified zone is much 
smaller than in the convective zone (figure |^) . We noticed a large difference 
in the stratified zone PSD at intermediate frequencies between the first 8 
hours of measurements (thin red line) and the following 64 hours (thick red 
line). But in both cases, no favorite frequency of wave propagation appears 
in the data analysis. We discuss the similar features of the two PSD curves 
below. 

Integrating the signals between frequencies 0 and N, we estimate the ratio 
between the wave energy and the convective energy to be about 1% at this 
depth, consistent with the estimates of Ansong and Sutherland (2010). The 
broad features of the stratified PSD follow those of the convective PSD, with 
most of the energy concentrated at low frequency and decreasing power with 
increasing frequency. However, there is a clear change of behavior around 
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Figure 3: (a) Normalized probability density function and (6) power spectral 
density (obtained with the Matlab function pwelch) of the temperature sig¬ 
nals shown in figure [^a), the data in black corresponding to the temperature 
fluctuations in the convective zone and the data in red corresponding to the 
temperature fluctuations in the stratified zone. In (6), the thin red line shows 
the power spectral density obtained from the first 8 hours of measurements, 
while the thick red line shows the power spectral density obtained from the 
following 64 hours of measurements. Typical power law dependences are 
indicated as dashed lines. The slopes corresponding to the convective sig¬ 
nal are estimated from the data published in Wu et ah (1990). The slopes 
corresponding to the stratified signal are indicative only. 
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the frequency 2 x 10“^ Hz. From the main convective excitation frequency at 
2 X 10“^ Hz to this frequency 2 x 10“^ Hz, the PSD slope in the stratified zone 
is significantly steeper than the slope in the convective zone. This is in quali¬ 


tative agreement with the analytical model of Lecoanet and Quataert (2013), 


which studied the waves produced by Reynolds stresses from Kolmogorov- 


type turbulence. As described by Lecoanet and Quataert (2013), the under¬ 


lying physical explanation is that large frequency excitation sources a priori 
correspond to smaller lengthscale fluctuation patterns: their e-folding depth 
of evanescence during their propagation in the convective zone is thus smaller. 
Additionally, one also expects smaller structures to be more rapidly damped 
by viscous and thermal diffusion, both in the convective and stratified zones. 
On the contrary, above 2 x 10“^ Hz, the PSD slope in the stratified zone 
is significantly less steep than the slope in the convective zone. We argue 
that this is directly related to the existence of propagating internal waves, 
and to their peculiar property that for a given horizontal wavelength, the 
damping length is an increasing function of frequency up to a; = 2N/ -^/S, as 


already noticed by Taylor and Sarkar (2007) in the case of waves generated 


by a turbulent bottom Ekman layer. 


As in Michaelian et al. (2002), we see a large difference of two orders of 


magnitude between the early and the late temperature PSD. We attribute 
this to more plume activity within the convection zone during the early times 
of the experiment, even if no sign of this randomness is seen in the statistics 
of the convective signals. Finally, the slope of the PSD of the first 8 hours 
of data (thin red line) steepens near the buoyancy frequency N. This is 
consistent with the transition between propagative internal waves expected 
for frequencies below A, and evanescent internal waves expected for frequen¬ 
cies above N. This is not observed in the PSD of the following 64 hours 
measurements, probably because of the limited precision of the temperature 
measurements. The spectrum then converges towards white noise, which can 
also been seen in the convective spectrum. 


4 Analysis of the velocity signals 


To perform velocity measurements in our set-up, the foam insulation is 
removed from the two long sides of the tank. As a result, the side heat 
losses significantly increase, and the typical depth of the convective zone is 
now 9 cm. Figure presents the main characteristics of the velocity field 
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(a) time t=437.5 s 



(b) log^Q(std of the kinetic energy) 
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Figure 4: (a) Example of the velocity field measured by PIV in the convective 
zone, with the color-plot on the left of the gray rectangle corresponding to 
three times the horizontal velocity and the color-plot on the right of the 
gray rectangle corresponding to the vertical velocity. Note that both signals 
come from the analysis of the same PIV movie taken on the right side of the 
tank. We Hipped the horizontal field across the central axis of the tank for 
visualization purposes. See also in supplementary material the corresponding 
movie “ConvectiveField.avi” showing the time evolution of those velocity 
measurements over 885 s. (6) Logarithm of the standard deviation of the 
velocity squared over 885 s (in (cm/s)^) and streamlines of the time-averaged 
velocity field, illustrating the mean organisation of the flow and the location 
of the main velocity fluctuations, (c) Example of the measured velocity 
fluctuations as a function of time: with reference to figure (a), the black line 
shows the vertical velocity at location x = 9 cm,y = 4.5 cm, and the gray 
line shows the horizontal velocity multiplied by —3 at location x = 5 cm, y = 
7.5 cm. 
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in the convective zone (see the corresponding movie “ConvectiveField.avi” 
in the supplementary material). As illustrated in figure |^a), the velocity 
field consists mostly of two symmetric large cells filling the whole depth of 
the convective layer, with strong, localized, rising cold plumes in the center 
and broader sinking motions on the sides. Figure |^6) shows the streamlines 
interpolated from the mean 2D fiow determined by averaging the velocity 
measurements for 885 s, starting from initial points regularly spaced in the 
x-direction on the horizontal lines y = 0.8 cm and y = 7.7 cm (i.e. the bottom 
and top of the PIV domain). This figure also shows the associated standard 
deviation of the kinetic energy, as a proxy for quantifying the velocity fiuc- 
tuations superimposed on the mean flow. Velocity fluctuations are mostly 
localized along the central axis of the tank, related to the passage of cold 
plumes forming along the lower boundary and then advected by the mean 
flow. The temporal evolution of the vertical and horizontal velocities at the 
locations where they are maximal is shown in figure |^c) . A typical period 
of about 75 s dominates both signals, corresponding to a frequency of about 
0.013 Hz. This main convective excitation frequency is significantly larger 
than the main convective excitation frequency of 2 x 10“^ Hz observed in 
the temperature measurement experiment in section This is probably due 
to the fact that the height of the convective layer is significantly smaller in 
the present case, while the injected heat flux is the same in both cases, corre¬ 
sponding to a temperature drop of about 4‘^C' through the lower conductive 
thermal boundary layer of about 0.6 cm (see Perrard et ah, 2013): we thus 


expect the turn-over timescale to be shorter for the smaller convective cells 
observed here. 


The temporal evolution of the wave field excited by this convective flow 
is measured by PIV during sequences of more than 10 minutes. A movie of 
measured wave velocities is provided in the supplementary materials (“In- 
ternalWaves.avi”). Figure]^ shows a frame of this movie, and figureshows 
the corresponding space-time diagram obtained by extracting data from the 
vertical line at x = 11.7 cm, so as to illustrate the wavy character of the 
velocity pattern. Figures [^a) andj^a) show the complete signal. Measured 
wave velocities are typically 25 times smaller than the convective velocities, 
and decrease rapidly away from the interface. The flow shown here corre¬ 
sponds to the superposition of waves traveling in all directions and with a 
wide range of frequencies. To better illustrate the wave patterns, we first 
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(a) full signal at time t=86 s 


X (cm) 


X 10 



(b) right-propagating signal filtered at f=0.013±0.001 Hz 



(c) right-propagating signal filtered at f=0.06±0.001 Hz 



Figure 5: (a) Example of the horizontal velocity field measured by PIV 
in the stratified zone, (b) and (c) show the corresponding signals includ¬ 
ing right-propagating waves only, respectively filtered around the frequency 
/ = 0.013 Hz (i.e., close to the typical convective frequency) and around the 
frequency / = 0.06 Hz (i.e., close to, but below, the buoyancy frequency). 
Note the changes in the color scale, 1-^ith interfacial velocities in panel (6) 
being more than one order of magnitude larger than interfacial velocities in 
panel (c). See also in the supplementary material the corresponding movie 
“InternalWaves.avi” showing the time evolution of these velocity measure¬ 
ments over 300 s. 






















(a) full signal at x=11.7 cm 
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Figure 6: Space time diagrams showing the time evolution along the vertical 
line at x = 11.7 cm of the horizontal velocity field measured by PIV in the 
stratified zone (see the corresponding movie “InternalWaves.avi” in the sup¬ 
plementary material). As in figure(o) shows the full signal, while (6) and 
(c) show the corresponding signals ijgluding right-propagating waves only, 
filtered around the frequency / = 0.013 Hz and / = 0.06 Hz, respectively. 











spatially filter the complete velocity signal to show only right-propagating 
waves. In our case the energy goes from the interface towards the top of the 
tank, so these correspond to wave vectors oriented towards the bottom right. 
We also filter our signal in time to only show waves in a narrow frequency 
band of ±0.001 Hz around a mean filtering frequency. In figures [^6) and 
1^6), we chose the filtering frequency / = 0.013 Hz, i.e. the typical frequency 
from the convective field (see e.g. figure]^). In figuresj^c) and|^c), we chose 
the filtering frequency / = 0.06 Hz, close to but below the mean buoyancy 
frequency of the stratified zone N ~ 0.07 Hz. 

The two types of waves previously described in the transient experiment of 


Michaelian et ah (2002) are clearly recovered in our data. Most of the energy 


is in waves with long horizontal wavelength and low frequency as shown in 
figures H*) and 1^6). Those waves concentrate close to the interface, and 
are rapidly damped as they propagate away from it. This region close to 
the interface corresponds to the buffer zone that was described by [Perrard 


et ah (2013), where the buoyancy frequency rapidly changes from zero in the 


convective zone to a constant positive value, and where low frequency waves 
propagate at angles as small as 10°. Additionally, some energy is carried by 
waves with shorter horizontal wavelengths and higher frequencies, as shown 
in figures [^c) and|^c). 


To further prove this last point, we show in figure spectrograms of the 
velocity signal as a function of height above the interface, with (a) showing 
the raw result and (6) showing the result normalized by the power spectral 
density at the interface. In both cases, a cut-off frequency close to the typical 
buoyancy frequency N ~ 0.07 Hz appears: energy is carried in the stratified 
zone by propagating waves only. Also, it clearly appears from the renormal¬ 
ized spectrogram in (6) that excited waves with frequency close to N are less 
rapidly damped, confirming our interpretation of the temperature spectrum 
in section]^ No favorite frequency of propagation appears here: the observed 
spectrum results from a competition between the way the injected wave en¬ 
ergy is distributed among the different horizontal wavelengths and temporal 
frequencies, and the way the waves are damped away from the interface. 


5 Conclusion 


We have re-investigated the “4°C experiment” first introduced by Townsend 


(1964), which allows a single self-organising system to produce a turbulent 
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Figure 7: Spectrograms of the horizontal velocity in the stratihed zone, mea¬ 
sured by PIV over 13 minutes (see an extract of the data in figures and 
1^ as well as in the movie “InternalWaves.avi” in the supplementary materi¬ 
als). Along each horizontal line in (a), the color plot shows the logarithm of 
the x-averaged power spectral density at a given depth, obtained using the 
Matlab function pwelch. (b) shows the same results when renormalized by 
the amplitude of the power spectral density at the interface (y = 0). 
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convective layer adjacent to a stably stratified one, as is encountered in atmo¬ 
spheric and stellar flows. We provide here a set of quantitative high-precision 
measurements that can help validate analytical and numerical models of such 
complex systems, whose dynamics range over wide time, length and velocity 
scales. The combination of spectral analysis of the temperature fluctuations 
and PIV velocity measurements allows us to describe the complex internal 
wave field excited in the stratified zone by the turbulence. In agreement with 


previous studies of closely-related configurations (Michaelian et ah, 2002; An- 


song and Sutherland, 2010), we describe the flow pattern as a superimposi¬ 


tion of low frequency, long horizontal wavelength waves close to the interface, 
and higher frequency, smaller horizontal wavelength waves propagating away 
from the interface. 

Previous studies also report the existence of a favorite angle of propa¬ 
gation, which was explained by Dohan and Sutherland (2003) as a resonant 
coupling between turbulent fluctuations and the excited internal waves that 
optimize the vertical transport of horizontal momentum. Such a selection 
is not observed in our experiment, where waves are of very small ampli¬ 
tudes: this conclusion supports a non-linear origin of the angle selection 
process in previous experiments, where larger amplitude waves were excited. 
In the present case, we rather hypothesize that the internal wave field is 
passively excited by the convective field, and that the observed wave field 
results from two contributing physical mechanisms: (1) the spatiotemporal 
characteristics of the turbulent excitation; and, (2) the selective damping of 
small wavelength waves and/or of low frequency waves during their propaga¬ 
tion in the stratified zone, as fully described by the dissipative, linear, plane 
wave theory (Taylor and Sarkar, 2007; Lecoanet et ah, 2015). The remaining 
question of the main process of wave excitation, i.e. interface deflections or 
bulk Reynolds stresses within the convective zone, is addressed in the pa¬ 
per Lecoanet et ah (2015), which studies direct numerical simulations and 
simplified numerical models of our experiment. Non-linear effects potentially 
responsible for zonal flow generation and for frequency selection will be ad¬ 
dressed in future work, involving a larger, three-dimensional, experimental 
set-up. 
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